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ABSTRACT.  In  the  field  of  spacecraft  attitude  control  and  stabilization, 
the  dual-spin  maneuver  has  an  important  place.  Here,  we  consider  a  problem 
which  is  a  multibody  analog  of  the  dual-spin  problem.  The  dynamical 
equations  are  derived  using  a  modified  form  of  the  Euler- Lagrange  equation 
on  the  special  orthogonal  group  50(3) .  It  is  then  shown  that  with  a  suitable 
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1  Introduction 


In  designing  a  communication  satellite  or  an  interplanetary  probe,  engineers  are  often 
faced  with  the  requirement  that  the  spacecraft  be  able  to  maintain  a  fixed  orientation 
relative  to  some  inertial  frame.  The  dual-spin  technique  is  a  simple,  commonly  used 
technique  for  meeting  this  requirement.  A  dual-spin  spacecraft  consists  of  the  spacecraft 
body  and  on-board  motor-driven  symmetric  rotors.  In  the  presence  of  a  suitable  damping 
mechanism  and  for  sufficiently  high  rotor  velocities,  the  attitude  acquisition  can  be 
achieved.  The  final  state  is  a  steady  spin  about  a  fixed  axis.  This  single  rigid  body 
dual-spin  problem  has  been  studied  extensively  before,  e.g.  [8]  [7]  [5]  and  the  references 
therein.  A  rigorous  proof  of  asymptotic  atability  can  be  found  in  [8]. 

In  this  paper,  we  study  a  multibody  analog  of  dual-spin  problems.  We  consider  a 
system  (see  Fig.  1)  consisting  of  two  rigid  bodies  connected  by  a  ball-in-socket  joint.  There 
are  also  three  rigid  symmetric  rotors  mounted  on  the  centers  of  mass  of  each  body.  One  set 
of  the  rotors,  called  driven  rotors ,  are  set  in  a  constant  relative  motion  to  the  carrier  body. 
The  other  set  of  rotors  carried  by  the  second  body  provide  damping  torques  to  the  overall 
motion.  These  are  called  damping  rotors.  We  shall  prove  that  with  an  additional  damping 
mechanism  on  the  ball-in-socket  joint,  the  motion  is  asymptotically  stable  and  approaches 
the  stable  equilibria  of  the  coupled  system.  In  deriving  the  dynamical  equations,  we  will 
use  explicit  representations  of  higher  order  tangent  bundles  of  the  special  orthogonal  group 
50(3). 

In  the  following,  we  shall  describe  the  system  configuration  and  compute  the 
Lagrangian  of  the  system  in  Section  2.  We  then  discuss  a  version  of  the  Euler- Lagrange 
equation  for  the  motion  on  the  special  orthogonal  group  50(3)  in  Section  3.  We  will 
generalize  this  approach  to  derive  the  dynamical  equations  for  our  problem.  In  Section  5 
we  shall  prove  an  asymptotic  stability  theorem  by  using  a  Lyapunov  function  and  invoking 
the  LaSalle  invariance  principle. 

The  problems  of  multibody  systems  in  space  are  very  complicated.  The  dynamical 
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behaviors  of  these  systems  are  highly  coupled,  highly  nonlinear,  and  poorly  understood. 
The  works  of  [13]  [10]  [6]  help  us  to  understand  the  phase  portraits  of  problems  such 
as  the  one  discussed  in  this  paper.  Also,  the  underlying  control  problems  associated  to 
multibody  systems  are  quite  complex.  Controllability  questions  for  rigid  bodies  with  and 
without  rotors  were  addressed  in  [5]  [3]  [4]  [2].  Similar  results  for  multibody  problems 
are  desirable.  Here  we  simply  set  up  the  Lagrangian  framework,  formulate  a  basic  control 
problem  and  then  prove  an  asymptotic  stability  theorem  for  a  multibody  system  under  a 
dual-spin  feedback  law.  The  questions  about  attitude  acquisition  and  rates  of  convergence 
are  under  investigation  and  we  hope  to  report  on  these  in  a  later  paper. 


2.  System  Configuration 

We  consider  a  system  of  two  rigid  bodies  connected  by  a  three-degree-of- freedom 
spherical  joint  with  rotors  mounted  on  each  body.  The  system  under  consideration  is 
shown  in  Fig.  1.  The  inertial  observer  is  at  the  center  of  mass  of  the  system  of  bodies. 

Ignoring  the  specific  kinematic  relationships  between  the  rotors  and  the  bodies,  the 
unconstrained  configuration  space  Qu  is  parametrized  by  the  attitudes  of  these  eight 
bodies, 

Qu  =  {  (Bi,  Si,  52,  S3,  B2,  D\,  D2,  D3)  }, 

=  50(3)  x  (50(3))3  x  50(3)  x  (50(3))3. 

To  account  for  the  body-rotor  relations,  we  have  the  following  constraints  between  the 
attitudes, 

Si  —  B\R(xi,6i),  i  =  1,2,3, 

Di  =  B2R(yi,  <f>i),  t  =  1,2,3. 

where  R(x{,9i)  is  the  rotation  about  the  x,  axis  by  the  angle  6, ,  e.g. 
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J 


Figure  1.  Two  Rigid  Bodies  with  Rotors 

(cos  83  sin  83  0\ 

—sin  83  cos  83  0  J  . 

0  0  l) 

With  these  constraints,  the  configuration  space  is  then 

Q  =  50(3)  x  (S1)3  x  50(3)  x  (51)3. 

Assume,  for  convenience,  the  centers  of  mass  of  rotors  5,  are  at  the  center  of  mass 
of  carrier  body  ,  and  the  centers  of  mass  of  rotors  Vi  axe  at  the  center  of  mass  of  the 
carrier  body  B2 .  It  can  be  easily  found  from  Fig.  1  that  we  have  the  following  kinematic 
constraints, 
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with 


rs.  =  n,  t  =  1,2,3, 

rD.  =  r2,  t'  =  1,2,3, 

r2  =  n  —  Bidi  •+■  i?2d2 
(mj  +  ms)  ri  =  -(m2  +  mD)  r2 

ms  =  mSl  +  mSj  +  rnS3, 
mD  =  mCl  +  mCj  mDs, 


where  mi ,  ms. ,  m2  ,  mDi  axe  the  masses  of  the  corresponding  bodies.  By  using  standard 
techniques,  the  total  kinetic  energy  of  the  system  can  be  written  as 

K  =  \  mil!  II2  +  \tr  {BJxBl)  +  i  m2||  f2  ||2  +  itr  (7?272£j) 

+  S  II2  +  ^MSJs.Sf)) 

+  ]C  ||2  +  ^trp.lD.Pn)  , 

where  7i ,  J2  ,  7s, ,  7#,  are  the  coefficients  of  inertia  of  the  corresponding  bodies.  Assuming 
there  is  no  potential  energy,  the  Lagrangian  on  the  tangent  bundle  to  the  unconstrained 
configuration  space  Qu  can  be  written  as 

L  =  \  tr  (B,  h  Bf)  +  \  !r(B,  J2  B\)  +  i  e  ||  B,d,  -  B2d2  f 

,  »  ,  .  .  .  .  .  (2-1) 

+  5  E  (tr(S.  Is‘  S‘~>  +  ,r(A  D]))  , 

»=1 

where 

_  (mi  -f  rn5)(m2  +  mD) 
mi  +  ms  +  m2  +  mD 


Now  we  let 

B\  =  BiQi,  B2  =  B2  Q2, 
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where  the  map  '  is  the  canonical  isomorphism  from  R3  to  so( 3) ,  the  skew-symmetric 
matrices.  Since  Si  =  BiR(xi,8i),  we  have 


S,  =  BiR(xi,8, )  +  BiR(xi,8, )  =  Bi(Cli  + 

and  thus 

ftsi  =  +  Si, 

where 

f? 

>  *2  =  J  &2 
\0 

We  can  then  write 

tr(S,Is,  ST)  =  tr({Cl1+st)R(xi,8l)ISiR(xt,0l)T(Cl1  +  st)T). 

We  naturally  assume  that  the  rotors  have  material  symmetry  about  the  axis  of  rotation, 
i.e. 

R(r„e,)Is>  R(xiA)T  =  Is,, 

and  we  get 

tr(S,  1st  Sf)  =  tr  +Si)ISi(& i  +5i)r)  =  (  +  s»>  +Si)  ). 

Similar  derivations  can  be  applied  to  the  rotors  T>i .  By  substituting  these  formulae  in 
(2.1),  the  Lagrangian,  L  :  TQ  — »  R,  can  then  be  written  as 
L(Bi,  0,,  B2,  <f>i,  fli,  0,,  ^2,  4>i,  *  =  1,2,3) 

=  —  (  ^1  )  +  —  (  fl2,  3%  ^2)  +  e  (  fli,  d\  B\  B2  d2  ^2)  (2.2) 

+  |(  0,7s©  )  +  («,,  7s©  )  +  |<  $,7D«j>  )  +  (  ft2,  7D*  ). 
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with 


3 

J,  =  I,  +  €%&+  ]TlS, 

t=l 

3 

h  =  h  +  t%4,+  £Id, 

i=1 

Is  =  diag((Is.)i  ,(Is,)2  ,(Is»)3), 
ID  =  diag((IDl),  ,(ID,)»  ,(Id,)s), 


where  (Is;  )j  denotes  the  j-th  diagonal  element  in  the  moment  of  inertia  matrix  Is, .  The 
physical  meaning  of  J i  is  that  it  is  the  total  moment  of  inertia  of  B\  plus  rotors  referred 
to  the  joint.  This  Lagrangian  will  be  used  later  as  a  basic  entity  to  derive  dynamical 
equations.  Our  approach  to  the  derivation  of  the  dynamics  is  novel.  We  illustrate  this  by 
a,  simplified  example  in  the  next  section. 


3.  Euler-Lagrange  Equation  for  Special  Orthogonal  Group  50(3) 

We  denote  an  element  in  the  special  orthogonal  group  50(3)  by  a  3  x  3  matrix  A 
satisfying  the  identity  AT A  =  I  and  det(A)  =  1 .  The  elements  in  the  tangent  bundle 
T50(3)  are  usually  expressed  as  (A,A£l).  This  is  a  global  representation  for  the  tangent 
bundle  of  50(3) .  On  the  other  hand,  we  recall  that  the  classical  Euler-Lagrange  equations 
are  based  on  local  coordinates.  In  [n],  Vershik  and  Faddeev  discuss  an  invariant  form  of 
the  Lagrange- D’Alembert  Principle  which  gives  the  classical  Euler-Lagrange  equation  with 
external  forces  in  local  coordinates.  Since  our  representation  for  T50(3)  is  global,  it  is 
useful  to  write  the  form  of  Lagrange- D’Alembert  Principle  in  that  representation.  We  now 
show  how  to  derive  the  dynamical  equations  for  motions  on  SO (3) . 

Lagrangian  Mechanics  is  about  second  order  equations.  Thus  we  need  to  find  a 
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representation  of  the  second  tangent  bundle  on  SO( 3).  We  first  make  the  following 
observation. 

d  -  •  -  *  * 

—t~aCi  —  Aft  -f*  Aft  =  A(flft  4“  ft). 
dt 

Thus  a  special  second  tangent  vector  must  be  of  the  form 

(A,Aft,Aft,A(ftfi  +  6))  . 

In  general,  let  (Au,  W)  €  A^TSO( 3).  In  order  to  have  this  vector  generate  a  curve 

(Ae^'Ae^ifl  +  tibj)  <E  TSO( 3), 

which  passes  through  (A,  Aft)  when  t  =  0,  we  must  have 

W  =  4-  Aet6(ft  +  fu>)  =  A(uCl  +  w). 
dt  (=o 

Any  element  in  TTSO( 3)  can  therefore  be  written  as 

^A,  Aft,  Au,  A(uft  +  u>)^  .  (3.1) 

Next  we  look  for  a  canonical  representation  for  an  element  in  T*TSO( 3).  First,  we  note 
that  the  trace  pairing 

<  A,  B  >  =  i tr(ATB ), 

provides  us  a  standard  way  to  define  the  elements  in  T*SO( 3),  i.e.  Ad  £  TASO{ 3),  and 

<  Aa,Au  >  =  j-tr(dT AT Au)  —  (  a,  u  ), 

where  (  a,  u  )  denotes  the  Euclidean  inner  product.  Letting  u /(A,  A)  =  (a,  ft)  € 

T('a,A,TS°( 3)' 

u>(A,  A)  ^Au,  A(uft  +  u>)^  =  ^tr(aT  Au)  +  ^<r(/?TA(uft  +  tf>)). 

In  order  to  have  a,  b  £  R3  ,  such  that 
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u>(A,  A)  (Au,  A(uCl  +  =  (  a,  u  )  +  (  b,  w  ), 


we  ask 

a  =  A(bCl  -j-  a) 

P  =  Ab. 

Thus  we  obtain  the  representation  for  any  element  in  T*TSO( 3), 

(A,ASl,A(bn  +  a),Aty  .  (3.2) 

These  global  representations  (3.1)-(3.2)  of  the  second  tangent  bundle  and  the  dual 
of  the  second  tangent  bundle  on  50(3)  prove  to  be  useful  in  finding  the  derivatives  or 
•variations  of  a  function  (Lagrangian)  on  T50(3)  and  in  deriving  the  reduced  Poisson 
bracket.  In  the  following,  we  state  the  Lagrange-D’ Alembert  Principle  in  terms  of  these 
representations.  Following  an  argument  similar  to  the  one  used  in  deriving  the  formula  in 
local  coordinates  (see  [1],  pp.  215),  we  get, 

THEOREM  3.1 

On  TSO( 3),  let  a  system  be  described  by  a  Lagrangian  L.  Lagrange-d’Alembert 
Principle  in  the  invariant  form  [n]  applied  to  the  motions  on  50(3)  gives  us  the  Euler- 
Lagrange  equation, 

<  D2L(A,a£i),  Au  >  =  <  DiL(A,Afl),  Au  >  +  <  a,  Aii  >, 

dt  (3.3) 

V  Au  €  7x50(3). 

where  a  is  the  external  force.  (Here  D\ ,  D2  are  the  usual  partial  Frechet  differentials.) 

I 

The  key  observation  here  is  that  with  these  representations,  the  space  of  vertical 
tangent  vectors  is  isomorphic  to  the  fibers  in  T50(3).  We  omit  the  details  of  the  above 
derivation.  The  only  difference  between  (3.3)  and  the  classical  Euler-Lagrange  equation  is 
in  the  explicit  interpretations  of  the  duality  pairing.  In  the  next  section,  we  will  use 
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the  procedure  of  this  section  as  a  model  for  working  out  the  dynamics.  We  remark 
that,  under  this  framework,  we  may  mix  the  local  coordinate  on  some  manifold  and  the 
parametrizations  for  50(3)  in  deriving  the  equations  of  motion.  This  is  very  helpful 
when  we  are  dealing  with  systems  which  have  a  mixed  configuration  space  (e.g.  Cartesian 
products  of  a  Lie  group  and  a  smooth  manifold). 


4.  Dynamical  Equations 

With  the  Lagrangian  formula  (2.2),  we  now  apply  Theorem  3.1  to  derive  the 
dynamical  equations  for  the  multibody  dual-spin  problem.  Let  L  be  the  Lagrangian 
function  expressed  in  terms  of  the  variables 

e  tq. 


First,  we  need  to  find  the  differential  of  L  in  a  form  analogous  to  (3.2).  It  can  be  found  by 
the  following  procedure.  Let  (t7j,  U2,  U3,  Ua,  Wi,  W2,  Wa,  WA)  G  T^Bi  e  Bi  *  ^  @  ^ ^TQ, 
which  can  be  written  as  the  form,  cf.  (3.1), 


(#1  Ui,  U2,  B2  U3,  4-tej),  U)2,  B2(u2&2  +  ^>3),  U'4). 


It  generates  a  curve  in  TQ  given  by 

(Bxe**1,  ©-Mu2,  B2e‘“2,  *  +  tu4,  Bie‘*l( fti+tuh),  0  +  tu>2,  B2e‘i3(fi2  +  fu>3),  $  +  *u>4). 


Here  we  have 
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dL-{ Uu  U2,  U3,  C/4,  Wu  W2,  W3 ,  W<) 


d_ 

dt 

d_ 

dt 


1=0 


L(B\etu l,  © -Mu2,  B2etUi,  $-ftti4, 

BietUl({l\  +  tw\),  &  +  iw2 ,  B2eiUi(Cl2  +  tw3),  $  +  <u>4) 
X(B1etttl,  ©  +  <u2,  B2etu*,  $  +  <u4, 

1=0 

flj -Mtiq,  0  +  <to2,  fi2  +  tu;3,  $-|-<u>4) 

OL  &L 

=  <  Db^iBiUi  >  +(  — ,  u2)+  <  -Db,.L,  j?2U3  >  +  (7^,  ^4) 

.  aZ  .  .dL  dL  .  .aZ  x 

+  W  Wl)  +  <5e’  W2)  +  W  “3)  +  {W  ^ 

The  canonical  form  for  dL(B\,  0,  B2 ,  $,  2?i,  ©,£2,  4>)  is,  cf.  (3.2), 

+  ai),  a2,  B2(b3Cl2  +  013),  a4,  B\bi,  62,  2?2Z>3,  • 

Let  A7i ,  iV2  be  given  by  the  formula 


we  have 


a  1 


fcj  = 


DBlL 

- 

dBsl  = 

B2^3, 

Nu 

dL 

03  =  ^3, 

az 

<*2 

a©’ 

a4  = 

a$’ 

dL 

1>2 

az 

l  dL 

az 

«v 

a©’ 

3  ao2 

,  64  = 

a$’ 

and  we  get  the  form  of  elements  in  (3.3), 


__  1  rw  a  91  _  ,  az  =  ~  dL 

DlL  =  0©’  B2{W2Q2  +  Nil  a¥ 

P2jL  1 5l0JV  a©’  S2dn2’  a4>  * ' 


(4.1) 


Next  we  model  the  external  force  a.  Let  TJ  be  the  torque  on  the  ball- in- socket  joint  in 
the  body  2  frame,  Ts, ,  T/5.  denote  the  torques  exerted  on  the  driven  rotors  and  damping 
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rotors  respectively.  Letting  B  =  B^B2,  the  shape  of  the  system,  we  can  write  a  in  the 
form 

a  =  (b^-BT3),  Ts,  B2T 'j,  Td)  . 


where 


/(TSl),\  /(TDl), 

(rs,)2  ,  TD  =  (TD,)2 

\(Ts,hJ  \(Td,)  3 


with  the  notation  (Ts,  )j  denoting  the  j-th  component  of  the  vector  Ts, . 

We  are  ready  to  apply  formula  (3.3).  The  first  component  is  worked  out  here.  The 
others  can  be  found  in  a  similar  way.  From 


B\ fl] 


dL 

dili 


+  B 


d  dL 
dt  dili 


-  B , 


dL 

dili 


Ch+Ni  +  (~BTJ)^j  =  0, 


we  get 


=  -«i  x  Mi  +  Ni  -  BTJ , 


where  Mi  =  .  This  can  be  rewritten  in  terms  of  Hi ,  etc.  Explicitly,  we  get  the 

dynamical  equations  of  the  system  in  terms  of  variables  in  TQ . 


Ji  iii 

-f-  0  tdiBd2il2 

=  -ill  X 

(Jlill 

+ 

ISQ  ) 

—  ediBil2d2H2 

-  BTJ , 

j2ii2 

+  +  td2BTdiili 

=  -il2  X 

(J2il2 

+ 

ID$  ) 

—  ed2BT  ilidiili 

+  TJ, 

is(ih  +  ©) 

=  Ts, 

iD(ii2  +  *) 

=  Td, 

Bi 

=  Biiih 

b2 

=  B7ii2. 

(4.2a) 
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For  the  realization  of  the  multibody  dual-spin  control  structure,  the  following  feedback  laws 


are  used, 


ts  =  Is  ii\. 

•*eT 

I 

II 

Q 

h 

(4.26) 

TJ  =  —  7(ft2-BTft1), 

where  /? ,  and  7  are  positive  definite  matrices.  The  first  equation  in  (4.2b)  makes  the 
relative  angular  velocities  between  the  driven  rotors  and  B\  be  constants.  The  other 
two  equations  signify  damping  torques  on  the  damping  rotors  and  the  joint.  With  these 
feedback  laws  and  the  following  transformation  of  coordinates  to  conjugate  momenta, 


Pi  =  Ji  Hi  +  t  di  B  d2  H21 

(4.3a) 

P2  =  Ise, 

(4.36) 

Ps  =  (h  -  1D)  ih  +  <  d,  BT  d,  fi,, 

(4.3c) 

P4  =  ID(Sl2  +  i), 

(4.3d) 

we  can  express  the  dynamical  equations  (4.2)  in  terms  of  pi  variables, 


Here  Oi ,  ,  and  4  are  the  expressions  of  ,  Sl2,  and  4  in  terms  of  p, ,  respectively. 

These  expressions  can  be  found  through  (4.3). 
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5.  Asymptotic  Stability 

Equations  (4.4)  describe  the  dynamical  behavior  of  the  system  under  investigation. 
We  are  now  ready  to  establish  an  asymptotic  stability  theorem.  First,  in  terms  of  the 
conjugate  momentum  variables  (4.3),  the  Lagrangian  can  be  written  as 

L  —  r>(  Pi  )  +  2^  ^S)_1P2’  P2  )  +  -(  ^2,  P3  ) 

+  2^  JP4’  P4  )  +  (  P2  )• 

We  define  the  function 

V  =  L  -  (flu  pi).  (5.1) 

It  can  be  shown  that  the  directional  derivative  of  V  along  a  trajectory  is 

dV  ~  ~  ~  ~ 

—  =  -(*,/?*)  -  (n2-BTShMn2-BTn1)). 

Thus  this  function  is  a  Lyapunov  function  candidate.  This  is  an  analog  of  the  core  energy 
used  in  [7]  to  justify  the  energy-sink  method  and  used  in  [8]  to  prove  asymptotic  stability. 
In  the  absence  of  damping  (  i.e.  7  =  /?  =  0  ),  the  system  is  a  hamiltonian  system  with 
hamiltonian  function  V" .  The  hamiltonian  structure  has  been  worked  out  in  the  ongoing 
Ph.D.  dissertation  of  L.-S.  Wang  [12].  It  is  a  generalization  of  the  one  derived  by  Grossman, 
Krishnaprasad,  and  Marsden  in  [6].  The  hamiltonian  structure  is  noncanonical. 

In  the  following,  we  shall  use  an  argument  similar  to  the  one  in  [8]  to  do  the  stability 
analysis.  It  can  be  shown  that  the  function 

||pi  +  P2  +  B(p3  +  P4)||2 

is  invariant  under  the  motion.  It  is  the  Casimir  invariant  of  the  hamiltonian  limit.  Since 
pa  is  a  constant  along  any  trajectory  (as  achieved  by  the  feedback  law  (4.2b)),  we  define 
the  momentum  variety 

=  {(pi,p3,p4,B)  €  R9  x  S0( 3)  :  IIpi  +P2  +  #(P3  +P4)li2  =  p2},  (5.3) 
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which  is  parametrized  by  p2  ,  and  p .  The  dynamical  motion  leaves  the  momentum  variety 
invariant.  This  will  be  the  domain  of  our  analysis  later  on.  Let 


K  ={  (Pi,P3,P.,B)€A/p'  :  —  =  0  }, 


(5.4) 


=  :  aDrV<  =  0,  ih-BTn,  =  0). 


r  2 


It  can  be  shown  that  7?.  is  parametrized  by  (p4,p,  2?)  as 

Pi  =  P  -  ( BJ2  +  tdx  B  d2)(7D)-1P4, 

p3  =  ed2  BT  d.  J-'p  +  (J2-Id  -e  d2  BT  dx  J-1BJ2)(JD)-1p4, 


(5.5) 


where  p  satisfies  ||p  +  P2||2  =  p2  ■  From  (4.4),  the  set  of  equilibria  in  M can  be  written 
as 

=  |(Pi»P3,P4,J5)  €  R9  X  50(3)  :  ||p  +  f>2||2  =  p2 , 

Pi  =  P  -  (B J2  4-  edx  B  d2){ID)~'p 4, 

p3  =  e  d2  BT  dxJxxp  +  ( J2-I°  -t  d2  BT  dxJ^BJ2){ID)-xp  4, 
(7d)’1P4  -  =  0,  Cl2  —  =  0, 

—  fix  x  (pi  +  p2  )  —  e  di  ilxBd2  Cl2  =  0, 

—  Cl2  x  (p3  +  p4  )  —  €  d2  Cl2BT  dx  =  0 

Now  we  characterize  the  maximal  invariant  set  in  1Z.  We  ask  that  the  directional  derivative 
of  the  functions  (7D)_1p4  —  fi2  and  Cl2  —BTQX  along  the  trajectory  be  0.  It  can  be  shown 
that  this  condition  is  equivalent  to 


fix  =  0,  f)2  —  0. 


From  (4.3),  we  have 

fpi\  _  (  Ji  ed.BdAfnA  _  (n A 
\p3J  Vc  d2  BT  dx  J2  -  ID  )  Vfi2  )  )  ' 


After  further  manipulations,  we  obtain  the  formula, 
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(—  Cli  x  (px  +  p2  )  —  e  d\  CliBd2  Q2  —  e  d\Bd2$l2 
—  CI2  x  (p3  -+-  P4  )  —  c  ^2  Cl2BT  d\  Q j  —  t  d.2BT d\Vl\  -(-  (3$ 

On  7?.,  B  =  0,  $  =  0,  we  thus  proved  that  the  maximal  invariant  set  is  exactly  the 
same  as  the  set  Y2  equilibria  for  (4.4).  From  LaSalle’s  theorem  [9],  we  know  that 

a  trajectory  of  motion  will  approach  the  maximal  invariant  set  in  the  limit.  We  can  thus 
conclude  the  following  theorem. 


±(&i 

dt  \&2 


THEOREM  5.1 

The  mechanical  system  (4.4)  asymptotically  approaches  one  of  the  stable  equilibria 
in  Y2  m,ps  5  or  equilibria  of  the  limiting  hamiltonian  system. 

II 

The  limiting  motions  of  the  system,  i.e.  the  equilibria,  are  exactly  the  same  notion 
las  the  relative  equilibria  discussed  in  [13].  The  techniques  for  finding  the  relative  equilibria 
there  can  be  applied  to  our  problem  by  additionally  including  one  gyroscopic  term  in  the 
augmented  potential  function. 


6.  Conclusions 

We  established  a  version  of  the  Euler- Lagrange  equation  for  the  motions  on  50(3) 
from  the  invariant  form  of  the  Lagrange- D’Alembert  Principle.  This  scheme  can  be  used  to 
derive  dynamical  equations  for  many-rigid-body  problems.  It  can  also  be  applied  to  certain 
continuum  mechanical  problems,  e.g.  the  special  Cosserat  rod.  We  applied  the  method 
to  a  multibody  analog  of  the  classical  dual-spin  problem  and  obtained  the  dynamical 
equations.  It  is  then  proved  that  with  a  damping  mechanism  on  the  ball-in-socket  joint 
and  the  damping  rotors,  we  have  asymptotic  stability.  The  driven  rotors  are  expected  to 
perform  the  role  of  attitude  acquisition.  This  problem  will  be  studied  further  in  the  future. 
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